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We investigate matter effects on highly relativistic neutrinos. The self-energy of 
neutrinos is determined in an electron or neutrino background taking into account 
resonance and finite width eff'ects of the gauge bosons. We find minor changes com- 
pared to the formerly used formula for the propagator function and large deviations 
of the effective width from the decay width of the gauge bosons considering higher 
moments of the electron or neutrino distribution function. 
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I. INTRODUCTION 

In the seminal work Wolfenstein explored the possibility of a change in the dynamics 
of neutrinos based on coherent forward-scattering with surrounding matter particles. Espe- 
cially interesting is the case of a flavor sensitive charged- current background leading to a 
change in the neutrino oscillation probabilities. The scattering amplitude was determined 
using the effective four-fermion Fermi interaction. In the picture of Wolfenstein, coherent 
forward-scattering leads to a refractive coefficient that accordingly changes the propagation 
of neutrinos in a medium similar to an effective potential of strength 

V = V2Gf{N,-N,), (1) 

where Gp denotes the Fermi coupling constant and N^. and A^g denote the electron and 
positron densities of the considered medium, respectively. 
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Of course, this picture is not valid for arbitrarily high neutrino energies. At an energy 
scale E ~ rriw ^10^ GeV, where mw is the W-boson mass, we expect the effective Fermi 
theory to break down and the dynamics of the W-bosons to become important. Energies of 
several GeV can indeed play a role in neutrino physics, if one, e.g., considers high energetic 
cosmic neutrinos [3] that could be accessible by experiments. However, this estimate is 
premature in the sense that the kinematics of the problem has been neglected. Arguing 
that the energy should be of order mw in the center-of-mass frame, one would come to the 
estimate (assuming the electrons to be at rest) 

s = 2Eu me + ml^^ 2E^ rrie ~ m^, (2) 

where rUe is the electron mass, which leads to 

E,^!^^ 10' GeV, (3) 
2me 

while for scattering with another neutrino (with mass of order 1 eV) 

1013 GeV. (4) 

Indeed, taking higher-order corrections in the neutrino energy into account, one finds to 
next-to-leading order for relativistic neutrinos P| 

8E, 

where (E^) and {Eg) denote the average energies of the electrons and positrons in the plasma, 
respectively. 

In this paper, we are interested in the limit of very high energetic neutrinos in matter, 
where the expansion used in Eq. (jS)) breaks down. In principle, one expects a resonance effect 
coming from the generation of real W-bosons by anti-neutrinos of center-of-mass energy 
s ~ The shape of the resonance depends on the finite decay width of the W-boson. 

This effect has been taken into account in Refs. by the educated guess that the finite 

width can be represented by just including an effective propagator of the W-boson into 
Eq. dH), leading to 

V = V2Gp [N, foi-sw) - N.foM] (6) 
with the propagator function given by 

fo{sw) = 1-. TTT 2~' (''') 



V = V2G1 



+ O {{El) , {El)) , (5) 
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where sw = '^E^rrie/m'^ and •jw = ^wl^^w denotes the width of the W-boson. 

The goal of this paper is to derive the propagator function from thermal field theory 
In addition, we give arguments how a resummation of the series in Eq. ^ changes our results. 
In the next section, we adopt the formalism of statistical quantum field theory as it was done 

riQ 

in Refs. p|, 121] and we present analytic as well as numerical results. Furthermore, we apply 
it to the cases with and without finite decay width. In a thermal background, a new energy 
scale enters (this scale is of order 1 keV if we consider neutrinos in the Sun up to order 
1 MeV for the cosmic neutrino background at decoupling) that will change the influence of 
the finite decay width of the gauge boson. Finally, in Sec. IIIH we summarize and conclude 
our results. 

II. THE NEUTRINO SELF-ENERGY IN A STATISTICAL BACKGROUND 

In the following, we carry out the calculations for the neutrino self-energy as was done in 
a temperature system in Refs. 0, 0, [l^ and applied to neutrinos in an electron background 
in Refs. . In contrast to these works, we expand in moments of the electron momentum 
(p") instead of the electron energy (-E" ) what is better suited for the investigation of the 
resonance regime. The corresponding Feynman diagram is depicted in Fig. Q The self- 
energy contributions coming from other species such as neutrinos or nuclear matter can be 
carried out along the same lines. 

W+{k-p) 




v{k) 

e"(p) 

FIG. 1: The Feynman diagram for the charged-current contribution to the neutrino self-energy. 



In the real-time formalism, the electron propagator reads 

1 



S{p) = {f + rrie 



— mj + le 

where F describes the electron background and is given by 



+ iT{p) 



(8) 



F(p) = 2tt6{p^ - ml)[e{p^)n+{p^) + ^( V)] (9) 
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and n±(p°) denote the particle distribution functions of electrons and positrons, respectively. 
The distribution functions are usually parameterized by the temperature and chemical po- 
tentials according to 

Furthermore, in the unitary gauge, the W-boson propagator is given by 

= , ~l ^. (9,. - ^1 • (11) 

The first term in the electron propagator is the usual zero temperature (T = 0) part, while 
the second term gives the following temperature contribution to the neutrino self-energy: 

X7^PiA^^(A;-p)|^7,Pi, (12) 

where the weak interaction coupling constant g is related to the Fermi coupling constant Gp 
as g'^/{Sm^r) = Gp/V^ and = (1 — 75)/2 denotes the left-handed projection operator. 
Performing the evaluation of the Dirac algebra, this can be transformed into 

^Ak) = --^ [{2ml, + k^HV'' + 2{k-V- mlQ)Jf\ Pp (13) 



with the definitions 



« - 1 (l^'"*^")"-'-"' «(-p>-(-/)i (,4;::^,, . (14) 

The expression can be split according to 

V^ = VuU^ + Vkk>', (16) 

where denotes the plasma four-vector that is (1, 0, 0, 0) in the rest frame of the plasma. 
The terms in the neutrino self-energy proportional to 1^ are related to the dispersion relation 
and the wave-function normalization and they give only corrections of order g'^ to the effective 
potential of the neutrinos ll|. In what follows, we only analyze the resonance behavior of 
the leading-order contribution to the neutrino self-energy which translates into the following 
effective potential for highly relativistic neutrinos 
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For simplicity, we neglect the mass (thermal and non-thermal) of the neutrino and we choose 
a frame such that k = {E^, 0, 0, E^). Out of the four integrations in Eq. |T^. three can be 
performed analytically, yielding 



Vik) 



3_ 

4 

+ 



dpp^ 



2ti^E,E, 



2n^EM 



n+{E,) 
n-{E, 



1 + 



w 



^pE, 



In 



ml, - ml + 2E^{Ee - p) 



1 + 



ml. - ml + 2E^{E^ + p) 
mlr - ml - 2E,{E, + p) 



ApE, 



V 



-ml- 2E^{Ee - p) 



(18) 



with Ee = \/ml + p"^. 

First, we will reproduce the result in Eq. (0). Thus, we expand the logarithm in the 
loop-momentum p 

mlr-ml + 2E,{E, T p) I , _ 4.E,p 



In 



ml,-ml + 2E,{E,±p) 



and obtain 

V{k) 

where we have defined 



9_ 
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as well as 



1 



-ml + 2Eyme 
dpp^ 



2n^ 



n+{Ee 



"e JO 

Note that in the limit 



-2E^m 
dpp"^ 



+ 0{{p)), 



27r2 



1 



e JO 



n-{Ee 



(p") = o, {E:) = m:, 



(19) 



(20) 



(21) 



(22) 



(23) 



this does not reproduce the result in Eq. (0), since there is an additional contribution of 
order {E~^) that can be neglected for high temperatures which was already pointed out in 
Ref. jsl. Including this contribution (and neglecting the positron background) yields 



AT, 



— ml 



1 



E,. 



{A{E^)-ml{E;')) 



(24) 



4 6(m^ - ml) 
which agrees with our result in the considered limit. 

So far, we have only expanded in the loop momentum p. The convergence of this expan- 
sion depends on the quotient of higher moments (p") and the resonance condition appearing 
in the denominators 



(m^ - ml - 2EymeY 



< 1. 



(25) 
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Hence, close to the resonance, our approximation breaks down and the divergence in Eq. ()2U|) 
is removed due to the finite values of the temperature T and the chemical potential fi in the 
electron distribution function. In Fig. |2l the full momentum dependence of the self-energy 
of neutrinos and anti-neutrinos is shown for the electron distribution in the Sun. 




FIG. 2: The propagator function f{sw) for neutrinos (left figure) and anti-neutrinos (right figure) 
without finite decay width. The electron background corresponds to the properties of the Sun. 



The height of the resonance can be estimated as follows. Suppose that the temperature 
T of the distribution function is not much larger than the chemical potential fi and notice 
that in this case the higher-order moments behave as (p") ~ /x". The potential as given in 
Eq. (j2(jp then depends on the center-of-mass energy s and the combination Ej^fi/m^. Being 
close to the resonance and identifying this parameter with the width as in Eq. ((7j), we obtain 
the estimate 

7^ = ^^^^r-V (26) 
ml 2m, ^30y' • ^^^^ 

The numerical value is 7^ ~ 0.0107, which is obtained by fitting the result with a resonance 
shape of the form in Eq. (jZj). Hence, the estimate in Eq. gives the right order of 

magnitude. 

Second, we would like to examine the inclusion of a finite width of the W-boson. The 
..a,.ve.e and longiu.dina, pan. of .he p.opaga.o.. ,„ E,. CD then change, to Q 

T ( ^ _ -1 ( 



KM = U.^ - ^ h (27) 

KM = V^^, (28) 
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where we have defined the W-boson width as 714/ = Vw /mw. Using the same arguments that 
led us to Eq. (fTTj) . we deduce that only the term proportional to the metric g^y contributes to 
leading order and that the longitudinal part of the propagator can be neglected. Analogously 
to Eq. (fTH|l . we obtain 



dpp^ 



a' 



2-K'^EyE^ 
dpp^ 



1 

Z AP Ey 



In 



ml^\1E,{E,-p)z 



m 



w 



1 

z Ap Ey z'^ 



In 



ml, + 2Ey{E^+p)z 
mlr-2Ey{Ee+p)z 



m 



w 



2E,{Ee-p)z 



(29) 



where we used the approximation m^r — ml ~ rri^ and the convention z = l+i'jw to shorten 
the expression. Performing similar steps as before and neglecting higher-order corrections 
in the temperature and the chemical potential of the electrons, we find 



V{k) = V2GF{Ne- Ne)^ 



m 



w 



m 



w 



E,,mp z 



(30) 



or equivalently, the W-boson propagator function is given by 

I - Sw 



(1 - Sw)^ -i 

in contrast to Eq. ((7j). Note that, in the limit sw 



(31) 



0, our result leads to the result of 
Wolfenstein, see Eq. (0), while Eq. (0) gives a correction by the factor 1/(1 + 7^). On 
the other hand, in the limit sw —>■ 00, we obtain an additional factor 1/(1 + 7^) in the 
effective mass of the neutrinos. This difference is due to the fact that we used a g^-dependent 
damping term in Eq. (j^Hj) . while a constant damping term can also be found in the literature. 
However, one has to have in mind that for the Standard Model value {jw — 0.0266) this 
correction is still out of the range of experimental testability, since absorption processes 
become strong close to the resonance On the other hand, the smallness of the width 
makes its influence compatible with the effect coming from flnite temperature and chemical 
potential, cf. Eq. (j^ . 

In Fig. El we have plotted the result of Eq. (jHTj) and the full numerical result including 
the effects of flnite temperature and chemical potential that can be obtained analogously to 
Eq. ()18|) . The numerical result corresponds to an effective width of 7efj ~ 0.0311. 

Thus, we arrive at the conclusion that in the Sun for the charged-current contributions 
the effect coming from the flnite width of the W-boson is larger than the temperature 



effect, i.e., 7eff 



7w 



> 



7^. This picture can change in different contexts. For example, 
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FIG. 3: The propagator function f{sw) with finite W-boson decay width. The plot fi{sw) 
shows the numerical result of Eq. H29|) including finite temperature and chemical potential effects, 
corresponding to 7efr ^ 0.0311, while the plot f2{sw) shows the result of Eq. (|3T|l . In both cases, 
the electron background meets the properties of the Sun. 

in supernovae, the rather large chemical potential of the electrons leads to 7^ ~ 1, and 
hence, the temperature effects dominate. Similarly, this effect is important considering the 
scattering with the neutrino background. Due to the small neutrino masses, the resonance 
appears for very high neutrino energies. Thus, using Eq. ()26p. we obtain 



2 E^^, 
7 = ~ 



(32) 



and the finite width from the Z-boson decay channels plays only a marginal role as long 
as the density of the neutrino background is larger than the mass scale of the neutrinos. 
This expectation is met in the numerical results of Fig. |3J where the resonance shape of the 
effective potential due to a neutrino background of temperature T = 10 eV is plotted for a 
neutrino of mass rriu = 1 eV and without chemical potentials /i± = 0. The resonance effect 
of the finite Z-width is almost completely washed-out, while the maximum of the resonance 
is shifted to smaller center-of-mass momentum. Since neutrinos are relativistic most of the 
time during the cosmological evolution (and one flavor could be relativistic even today), we 
expect the change of the effective width to be relevant for neutrinos from cosmic sources 
close to the resonance. 
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FIG. 4: The propagator function f{sz) with finite Z-boson decay width. The plot f{sz) shows 
the numerical result of Eq. H29|) including finite temperature effects (T = 10 eV) of the cosmic 
neutrino background (CNB) of a neutrino of mass rui, = 1 eV without chemical potentials fi± = 0. 

III. SUMMARY AND CONCLUSION 



We have analyzed the rare scenarios, where the neutrino energy is large enough such 
that the energy and momentum dependence of matter effects could be relevant. This is, 
e.g., the case for neutrinos from cosmic sources that has been analyzed in Ref. In 
this limit, we have derived the shape of the resonance from thermal field theory and found 
a small deviation from the formerly proposed propagator function. In addition, we have 
examined the change in the effective width due to the finite temperature and chemical 
potential and found that the effective width will be much larger than the decay width of the 
W-boson in supernovae or the Z-boson width in the case of scattering with the relic neutrino 
background. Finally, we note that the Mikheyev-Smirnov- Wolfenstein (MSW) Q, Q, Q 
resonance condition for effective mixing of neutrinos in matter will, in principle, be modified 
due to the temperature effects, since these effects can be interpreted as flavor conserving 
non-standard Hamiltonian effects 16]. However, the change will be completely negligible 
for solar and supernova neutrinos. 
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